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Setti and Woltjer have recently discussed a universal model for 

[2 1 

cosmic rays. This model i ased on an earlier paper in which the 
same authors assumed that the infrared background at submillimeter wave- 
lengths^ is universal, resulting from the superposed contributions 

of the infrared luminosities of Seyfert galaxies at large redshifts. 

If the cosmic ray and infrared outputs of these galaxies are about equal, 
the cosmic ray energy density in metagalactic space is comparable to the 
energy density in infrared photons. Since the latter could be as large 


as a few (eV cm*^) 


16,7] 


, this model, in principle, is capable of account- 


ing for the observed energy density of cosmic rays near earth. 

One of the principal observable features of the cosmic radiation is 
its energy spectrum. Since cosmological effects leave the cosmic ray 
spectrum invariant at relativistic energies only, it is of some interest 
to investigate -he effects of an expanding universe with evolving sources 
on low energy cosmic rays. 

Let (_ 0 /t) m q(p) dp be the volume emissivity of cosmic rays 
measured in particles p.r cm per second at some epoch t, where t Q is the 
present epoch and p is the particle momentum at epoch t. The function 
(t 0 /t) represents possible evolutionary effects, with m ranging 


from 0 (no evolution) up to perhaps m = 6.5 


[ 2 ] 


The expansion of the universe leads to a variation of p with time 

[ 8 ] 


which is given by 


o(t)R(t) * C.O'ns't 


( 1 ) 


where R (t) is the time-dependent scale factor of the universe determined 
from Einstein's field equations. If cosmic rays are freely moving particles 
in metagalactic space (i.e., the effects of matter and magnetic fields are 
negligible), the differential cosmic ray intensity, I(p 0 )» resulting from 
particle production in the space-time element dt dV can be written as 

tftto) q* 




where p * p c R (t 0 )/R(t), ^ * r/R is a dimensionless radial coordinate, 

dV * R 2 (t)7] 2 R(t)dudS, and 





( 3 ) 


is the invariant distance between two points embedded in the metric and 
traversed by a particle of velocity . Since u does not depend 

on the expansion of the universe, we have that 
Furthermore, , so that equation (2), inte- 

grated over all |plci u time, becomes 


MJU1IW V< U UUWO ftIV W 

it /dt 0 -p(i.ypw R(t)/«(^ 



( 4 ) 
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where and "t m is the epoch at which cosmic-ray sources began 

to be formed. 

Except for the factor , equation (4) is the same as the corres- 

* [9] 

ponding expression for photons . In particular, if \(p w p > 

I (po) oc po Po" n , i.e., a power law spectrum is modified by the 

factor at low energies. This result is independent of cosmological 

model or evolutionary effects. 

An interesting departure from this simple behavior could be caused 
by spatial inhomogeneities in the cosmic ray source distribution. Let us 
assume that there is a minimum invariant distance u min , such that at all 
epochs t, t m < t < t 0 , there have been no sources between us and u m; ^ n . 

In this case, the integral i- equation (4) has to be cut off at an upper 
limit t r < t 0 . The time t r is momentum dependent and is determined by 
solving the integral 


Um . , f 

J m •> 

t T 


(5) 


In order to associate this cutoff with an observable quantity, we consider 
photons emitted by a source at 1 %^ with redshift z m £ n . 


Um -» ~ f o <f+i) • 


( 6 ) 


The quantities dt/dz in equation (6) and R(t) in equations (4) and (5) 
depend on the cosmological model. In a zero-pressure Friedmann universe 

Jt/Ji * - if Ho ? 





•IfiilHr.ln'Tilafi-l ■ vn ..-i. 
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where n 0 is the mean matter density, and H 0 is 

Hubble's constant at the present epoch. Furthermore, as discussed by 

Setti and Woltjer^, a universal theory of cosmic rays 4a tenable only 

if the metagalactic density is less than about 10"^ cm”^, since otherwise 

the diffuse gamma ray background at £ 100 MeV would exceed the observed 

Clll 

upper limits . We have to use, therefore, a low-density model for 
which O and * * 1 /MoQtlZ) ^ . By changing the variable 

of integration in equation (4) to h * R(t>)/m - / and by combining 
equations (5) and (6) we find that 




m-/ 


where 


tc 


(f. + /?<>/ J (‘+ 2 ™*) 2 ! _ 4 

Zf>o ( fe + tyo +i ')((+ *•**) 


(7) 


( 8 ) 


The quantity z r is given in Figure 1 for various values of p G and 
*min- ^ f>«/ , if * 2 mm ( / + , whereas 

for f»f f . The value of the cutoff momentum p c is obtained 

by equating z r and z max . In the non-relativistic region 



/ /f St h tl'w 
( /+ 


( 9 ) 


Al.tP.'fjgtBL-’.dl 
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, 2 , mm4L * 2 




We have evaluated equation (7) for 
and * 0, 0.5 and 1. The results are shown in Figures 2 and 3, for 

m = 0 and 6.5, respectively. As can be seen, the low energy cutoffs are 
independent of m, but above the cutoff, the spectrum steepens with increas- 
ing m because, by increasing the evolutionary effects, the number of 
relativistic particles that are redshifted to low energies is correspond- 
ingly increased. 

[ 2 ] 

Setti and Woltjer have evaluated the energy density UJlf of 
the far infrared background due to Seyfert galaxies in a cosmological model 
with <^arf/2,. We can get similar results for ^.p* ° by using equation 
(7) with J3 0 «r/ and , where JO* and 

are the assumed density and luminosity of Seyfert galaxies at the present 
epoch. The energy density in infrared photons is then given by 




vn-3 


( 10 ) 


We have evaluated J° 0 = /o‘ 77 6W J ; L 0 ~ZS*lo\ r # 

and z max 58 2.5, and various values of m and z m £ n » The results are given 
in Table 2. As can be seen, even though the variation of with z m i n 

is large for m * 0, it becomes negligible if the evolutionary effects are 
large. The same effect can also be seen by comparing the spectra in 
Figures 2 and 3 at relativistic energies. Since on energy grounds, a 
universal theory for the infrared background as well as for the cosmic 
rays requires large evolutionary effects, we conclude that independent 
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of the energy densities in both infrared photons and cosmic rays, a 
spatial inhomogeneity in the source distribution will produce a low 
energy cutoff in the spectrum of the cosmic rays. 

Cosmic ray spectra with low energy cutoffs are not inconsistent 


with direct observations near earth. Goldstein et al. 


[ 12 ] 


have shown 


that because of energy loss in the interplanetary medium, cosmic ray 
particles below ~ 100 Me V/ nucleon are prevented from reaching earth. It 
is, therefore, not possible to sample the cosmic ray spectrum at low 
energies near the orbit of earth, and the cutoffs derived in this paper 
may in fact exist. 

Although there are no compelling arguments for a metagalactic model 
of cosmic rays, such a model cannot be ruled out on the basis of existing 
observations. A possible exception is the issue of the ionization and 
thermal state of the interstellar medium. In equilibrium, the rate of 


ionization of the interstellar gas ^ diould net exceed about 2.5 x 10 


-15 


(H atom sec) 


_][ 13] 


. For a cosmic ray intensity with the spectral shape 


given in Figures 2 or 3 for z m £ n * 0 and normalized to direct measurements 

at earth at high energies, ^ 5 x 10“ (H atom sec) \ where g is 

kinetic energy in MeV. In order not to overheat the interstellar medium, 

a cutoff is required at about 0.2 MeV. If this cutoff is produced by the 

mechanism suggested in this letter, from equation (9) with z max 85 2.5, 

we find that z min ■ 0.05. A higher value of z m £ n would be required if 

-2 5 

the source spectrum is steeper than p . Conversely, if the source 
spectrum is flatter or possesses an intrinsic low energy cutoff, z m £ n 
could be smaller than 0.05 or could in fact be zero. 


7 


Finally, if the bulk of the cosmic ray are metagalactic but are cut 

off at low energies, galactic supernovae could still produce large fluxes 

of low energy nuclei which would heat and ionize the interstellar medium, 

[14] 

as discussed recently by Ramaty et al. . In addition, in a metagalactic 

model of cosmic rays, the role of galactic supernovae would be limited to 

the production of cosmic electrons (and nuclei to about 1% of their 

observed energy density) and possibly, to the generation of certain short 

[15] 

lived ultraheavy nuclei which would be present in the cosmic rays 


We are grateful to Dr. E. A. Boldt for valuable discussions on the 


cosmological aspects of this problem. 
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The variation of 2 r , defined in equation 8, with momentum 
Differential energy spectra for m * 0. 

Differential energy spectra for m * 6.5. 
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FIGURE 3 


